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Abstract 

We investigate the influence of the universal extra dimension on the branching ratio 
in the B TT{p)i^i~ decay. Taking l/i? ~ {200— lOOOjGeV with one universal extra 
spatial dimension, which is consistent with the experimental data for B{B — > Xg^), 
B{B — > K* and the electroweak precision tests, we obtain that for both (/x, r) 

channels the branching ratio strongly depends on the compactification radius \/R. 
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1 Introduction 



The standard model(SM) has been successful theory in re-producing almost all experimental 
data about the interaction of gauge bosons and fermions. However, the SM is not regarded 
as a full theory, since it cannot address some issues i.e. , gauge and fermion mass hierarchy, 
matter- antimatter asymmetry, number of generations and so on. For these reasons, the 
SM can be considered as an effective theory of some fundamental theory at low energy. 

Extra dimension model [1] is one of the candidates trying to shed light on some of those 
issues. It can be categorized in terms of the mechanism of new physics where the SM fields 
are constrained to move in the usual three spatial dimensions(D3 bran) or can propagate 
in the extra dimensions( the bulk). The last one can be categorized as non-universal 
extra dimension(NUED) and universal extra dimensions(UED). In the non universal model 
the gauge bosons propagate into the bulk, but the fermions are confined to bran. In 
contrast, the UED allows fields to propagate into the bulk. The UED can be considered 
as a generalization of the usual SM to a D^j^n bran where N in the number of the extra 
dimensions [3]. The model introduced by Appelquist, Cheng and Dobrescu (ACD) [2] is the 
most simple example of the UED where just single universal extra dimension is considered. 
This model has only one free parameter in addition to the SM parameters and that is the 
compactification scale R. Mass of the Kaluza-Klein(KK) particles are inversely proportional 
to R, then, at large value oi 1/R the SM results can be almost recovered, since the KK 
modes, being more and more massive, are decoupled from the low-energy SM. 

Two types of study can be conducted to explore extra dimensions. In the direct 
search, the center of mass energy of colling particles must be increased to produce Kaluza- 
Klein(KK) excitation states, where KK excitation states are supposed to produce in pair by 
KK number conservation. On the other hand, we can investigate UED effects, indirectly. 
The indirect search at tree-level, where KK excitations can contribute as a mediator, is 
suppressed by KK number conservation. On the contrary, the same states can contribute 
to the quantum loop level where the KK number conservation is broken. As a result, flavor 
changing neutral current(FCNC) transition induced by quantum loop level can be consid- 
ered as a good tool for studying KK effects. The collider signatures and phenomenology of 
UED have been studied by Ref. [3] and [4, 5], respectively. These studies have provided a 
theoretical framework to examine some inclusive and exclusive decays with the ACD. 

FCNC and CP-violating are indeed the most sensitive probes of New Physics (NP) 
contributions to penguin operators. Rare decays, induced by flavor changing neutral current 
(FCNC) h s{d) transitions is at the forefront of our quest to understand flavor and the 
origins of CP violation asymmetry (CPV), offering one of the best probes for NP beyond the 
Standard Model, in particular to probe extra dimension. In this regard, the semileptonic 
and pureleptonic B decays have been studied with UED scenario [4] -[9]. They have obtained 
that the inclusive and exclusive semileptonic and pureleptonic decays are sensitive to the 
new parameter coming out of the one universal extra dimensions i.e., compactification scale 
1/R. 

New physics effects manifest themselves in rare decays in different ways: NP can con- 
tribute through the new Wilson coefficients or the new operator structure in the effective 
Hamiltonian, which is absent in the SM. Also, it may modify the SM parameters such as 
masses and CKM matrix elements. A crucial problem in the new physics search within 



flavour physics in the exclusive decays is the optimal separation of new physics efi^ects from 
uncertainties. It is well known that inclusive decay modes are dominated by partonic con- 
tributions; non-perturbative corrections are in general rather smaller[10]. Also, ratios of 
exclusive decay modes such as asymmetries for B ^ K{ K*, p, 7) decay [11]-[15] are 
well studied for new-physics search. Here, large parts of the hadronic uncertainties par- 
tially cancel out. The universal extra dimension with only one universal extra dimension 
belongs to the classes of NP, where the Wilson coefficients is modified by KK contributions 
[4, 5] in the penguin and box diagrams. Also, CKM matrix elements and masses are af- 
fected by ACD model. Obviously these modifications will affect the physical observables. 
In this connection, we try to investigate the effects of one universal extra dimension on the 
branching ratio(Br) of the B — > T^{p)t^£~ decay. 

The paper encompasses three sections: In section 2, we recall the effects of the UED 
on the inclusive h — > di^£~ decay and the general expressions for the matrix element and 
branching ratios of 5 — > 7r{p)i'^£~ decay is presented. In section 3, we investigate the 
sensitivity of Br and CP asymmetry to the copmactification scale(i?) and conclusion. 

2 Matrix element b dtl' in ACD model 

The QCD corrected effective Lagrangian for the decays h d£^i~ can be achieved by 
integrating out the heavy quarks and the heavy electroweak bosons in th SM: 



M{b 



s£H-) = ^V,,V. 



td 



■'9 
+C10 



diai^^^Rb 
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where are Wilson co-effeicents calculated in naive dimensional regularization (NDR) 
scheme at the leading order (LO), next to leading order (NLO) and next-to-next leading 
order (NNLO) in the SM[16]-[22]. 

As we mentioned above, in ACD model the new physics comes through the modification 
of the Wilson coefficients and the operator structures remain the same as SM. Considering 
the KK modes effects in the penguin and box diagrams, the above coefficients have obtained 
at NLO [4, 5]. Clearly, they depend on the additional ACD parameter i.e., R. These 
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coefficients can be expressed in terms of the functions F {xt,l/ R), Xt — -^i-, which is the 
generahzation of the corresponding SM function Fq (xt) according to: 



F {xt, 1/R) = Fo (xt) + J2Fn {Xt, Xn) 



(2) 



n=l 



with Xn = 

E{xt,l/R) 



M. 



iTLn = ^ [2] corresponding functions are C{xt,l/R), D(xt,l/R), 

1/R) and E' {xt,l/R), respectively. The Wilson coefficients in terms 
of these functions computed in [16]-[22]. We recall the formulae for the Wilson coefficients 
where we use the Wilson coefficients at the renormalization scale /x = m;, = 4.7 GeV. 
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w 

D' {xt, 
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where 77 = ^ilMnl^ and 



C,^^^(/.^) = 1, = —D'ixt, -^), C3^^^(/x.) = -^^'(x,, i); (4) 

where the wilson coefficients have been calculated in the leading order approximation. 
Moreover: 

14 16 6 12 

a-t = — tto = — a-i = — ad = 

23 23 23 23 

as = 0.4086 cte = -0.4230 ay = -0.8994 ag = -0.1456 
/ii = 2.996 h2 = -1.0880 /ig = /i4 = 

/i5 = -0.649 he = -0.0380 /i7 = -0.0185 hs = -0.0057. (5) 

The functions D' and E' in Eq. (5) are given as: 

. (8x^ + 5x^-70;,) a;g(2-3x,) 
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D'n{xt,Xn) 
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Xfii^ 7xn + 3Xyj) x^ , , 

+ 6 ^"TT^ 
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12(xi - 1) 

+ -^x„(l + x„)(-l + 3x„)ln-^^ (9) 

(1 + Xn){Xt + 3x2 _^ ^2^3 _ ^^(^^ ^_]^Q _^ Xn + X^ 
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With regard ACD model and in the NDR scheme we have 

CM = pr"" + ^^fei^ - 4z(x„ i) + p^£;(x„ i) (10) 

sm K K 



where P(f^^ = 2.60 ± 0.25 [20] and the last term is numerically neghgible(PE ~ lO^^). 
Besides 

^ n=l 
-j^ 00 

^(^t'o) = Zo(xt) + y2 Cn(xt, Xn) (11) 

n=l 



with 



Zo{xt) = 



18x^ - 163xf + 259j;? - lOSx* 



144(xt - 1)3 
38xf + 15x^2 
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+ [^M^^M±M^_l]ln., (12) 



Cn(a;t, a;„) = ^^^"T^p -8xt + 7 +{3 + 3xt + 7x„ - XtXn) In y^"^] (13) 

On the other hand, the effective coefficient is scheme independent. It can be parametrised 
as follows: 



where 



and 



a^^ = ei + At„6, (14) 

= Cg + 0.138a;(s) + g{rhc, s)Co{7hb) - ]^g{rhd, s)(C3 + C4) 

- ^^(m;,, s)(4C3 + 4C4 + 3C5 + Ce) + ^(BCg + C4 + 3C5 + Cg) , 
6 = [9{rhc,s) - g{mu,s)\{3Cx + C2) , (15) 
where rhq = mq/rrib, s = g^, Co(/i) = 3Ci + C2 + 3C3 + C4 + 3C5 + Cq, and 

2 4 2 5 -I- 4s 

2s(l + s)(l-2s) 5 + 9s-6s2 

~ 3(l-s)2(l + 2s) + 3(l-s)(l + 2s) ' ^ ^ 

represents the 0{as) correction coming from one gluon exchange in the matrix element 
of the operator Og [19], while the function g{mq,s) represents one-loop corrections to the 
four-quark operators Oi-Oq [21], whose form is 

8, X 8 4 2,^ 
g(mq, s) = -- ln(mj + ^ + g^? " 9 (2 + Vq) 



^\l-yJe(l-yq) 



+ 9{yq - 1) arctan | (17) 



where Uq = Arh'^/s. It should also be cmaphasized that Cf^ is in particular sensitive to the 
rhc in the NLO. To reduce this dependency NNLO calculation is necessary. 

Although long-distance effects of cc bound states could contribute to Cg^, for simplicity, 
they are not included in the present study. 



•Cio 

Cio is II independent and is given as: 



We aim to calculate the decay rate for B 7r{p)i~^i~ using the experimental allowed 
region for 1/R from the B{B — > Xgj) and B{B — > K*ix^ii~) decays. One has to sandwich 
the inclusive effective Hamiltonian between initial hadron state B{pb) and final hadron 
state 7i{pTr){p{pp)) to obtain the matrix element for the exclusive decay B 7r(p) 
Following from Eq. (1), in order to calculate the decay width and other physical observable 
of the exclusive B —>■ TT{pT:){p{pp))i'^i~ decay, we need to parametrize the matrix elements 
in terms of formf actors. 



2.1 Decay rate for B ^ ir i'^i decay 

The exclusive B 7ri~^£~ decay, which is described in terms of the matrix elements of the 
quark operators given in Eq. (1) over meson states, can be parametrized in terms of form 
factors (/+/~and fv). 

{7r{p,)\d^,{l-^')b\B{pB)) = r{q'){p.+PB)^ + f-{q')q„ (19) 



{7:ip,)\dia^,q-'{l + j')b\BipB)) = [q^ip. + Pb)^. - qMl ' ml)]Mq'), (20) 
Now, we can obtain the matrix element as: 



M^-^'^ = |^l-,,l-,^|(2Ap^^ + S0^7/ + (2G'p^^ + L'0%7'n, (21) 



where 



A = Cff+-2mBCff,, (22) 
B = Cf{r + n + 2'^CfMml-ml-q% 

G = C.of^, 

D = Cio(/+ + r), 

Prom this expression of the matrix element, for the unpolarized differential decay width 
we get the following result: 



= lmlK{3-v^){\A\^ +\G\^) + 16m^^r^\G\^ + irr4s\D\^ (24) 
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+ 8mj{l-r^- s)Re[GD*], 



with = ml/m%, A^r = + (s — 1)^ — 2r^(s + l),v = y^l — ^ and i = mi/niB- We use 
the results of the constituent quark model [23], where the form factors /t and /+ can be 
parametrized as: 

(l-gVT|)[l-aigVM2 + a2gV^'] ■ 
In this model, /_ is defined slightly different and it is as: 

^^"^^^ " [1 - (7igVM2 + ^2^7^^] ■ ^^^^ 
The parameters /(O), can be found in Table 1. 





/(O) 


0-1 


0-2 


/+ 


0.29 


0.48 






0.29 


0.76 


0.28 


/. 


0.28 


0.48 





Table 1: i? — > tt transition form factors in the constituent quark model. 



2.2 Decay rate for the B ^ p t^i decay 

Similar to the B ^ tt decay the following matrix elements defined in terms of form- 
factors must be computed for the B ^ p i'^£~ decay: 

{p(p„e)\d^,(l-^')b\B(pB)) = -6^.A.£>X ^^i^'^ -ie;(mB + m,)Ai(Q') 

TYIB + Tflp 



rriB + rup 

+ tq,{e*q)^[A,{q') - A,{q% (27) 

{piPp, e)\dia^,q''{l ± ^')b\B{pB)) = 4e^,^,e^*pyT,{q^) ± 2i[e;{ml - mj) 

- ipB+Pp),{e*q)]T2{q') 

± 2z(£*g)L-(pB+l>p)M^^^)T3(g'), (28) 

(p(p„£)|J(l + 7^)6|S(p^)) = -^^e*q)Aoiq% (29) 



where Pp and e denote the four momentum and polarization vector of the p meson, respec- 
tively. 

From Eqs. (27,28,29), we get the following expression for the matrix element of the 
B pe+i-decay: 



+ 



piq^D + e;E+{e*.q){pB)F 
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(30) 



(31) 



From this expression of the matrix element, we get the following result for the differential 
decay width: 

(dVP\ GW , ,,2 , A—, 



'■pi 
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4m|s|A|2 + ^(l + 12^)|5|2 



rriBrp 



^SX 1^12 2 



+ -^X^\G\' + -{l-rp + s)Re{B*G) 



+ imlXpis - 4t^)\D\' 
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2sr/^ 
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{2r + s){Xp + 2s + 2rp) - 2{2t\rp + 5s) + s{rp + s)} 



s , I 2Ap 



+ 3-A,|G|^ + 



sr. 



- 2i\rp - 5sO + (2^2 + s) - s(r^ + s) 



Rc{E*F) 



+ -^^Xp Re{G*E) + ^^A^(l -rp + s) Re{G*F) 

rUBTp Tp 



(33) 



with = mym\, Ap = + (s — 1)^ — 2rp(s + 1),'^ = ^Jl — ^ and t — me/m.B- The 
definition of the form factors are(see[24]): 

V{0) 



^(g^) = Ai(0)(l- 0.023^2), 
^2(?') = >l2(0)(l + 0.03V), 

^3(0) 



^^^^^ - 1-^75.3^' 
T2(g') = r2(0)(l-0.02g2), 

= T3(0)(l + 0.005g'). (34) 

with V{Q) = 0.47, ^i(O) = 0.37,^2(0) = 0.4,Ti(0) = 0.19, r2(0) = 0.19, T3(0) = -0.7. 



3 Numerical analysis 

In this section, we study the dependence of the total branching ratio on the compactification 
parameters (1/i?). The main input parameters in the calculations are the form factors. We 
use the results of Refs. [23] and [24] for i? ^ tt and S — > p transitions, respectively. Also, 
we use the SM parameters shown in table 1: 



Parameter 


Value 


ois{mz) 


0.119 




1/129 


mw 


80.41 (GeV) 


mz 


991.18 (GeV) 


si'n?{9w) 


0.223 


mi, 


4.7 (GeV) 




0.106 (GeV) 


m-r 


1.780 (GeV) 



Table 2: The values of the input parameters used in the numerical calculations. 

The allowed range in the ACD model for the Wolfenstein parameters at 1/R — 200GeV 
are: 0.076 < p < 0.260 and 0.305 < fj < 0.411 [4]. Note that, there is a small discrepancy 
with respect to the SM values in the higher values of 1/R. In the present analysis, they are 
set as p = 0.25 and fj = 0.34. In the Wolfenstein parametrization of the CKM matrix, Xtu 



is: 



p(l — p) — rf — ir] 
(1 - p)2 + 772 



+ 0(A2). 



(35) 



Furthermore, we use the relation 



(36) 



where A = sin^^vi^ • 

From exphcit expressions of the physical observables, one can easily see that they depend 
on both s and the copmactification radius (1/i?). One may eliminate the dependence on 
one of the variables. We eliminate the variable s by performing integration over s in the 
allowed region. The total branching ratio is defined as 



The branching ratio given by Eq. (37) depends on compactification radius(R). The conser- 
vation of KK parity (—1)-' , with j the KK number, implies the absence of tree level contribu- 
tion of KK sates at low energy regime. This allows us to establish a bound: > 250GeV 
by the analysis of Tevtron run I data[2]. The same bound can be obtained by the anal- 
ysis of measured branching ratio of 5 — > Xg^ decay [4, 5]. A sharper constrain on 
can be established by studying the zero point position of forward-backward asymmetry of 
B K*£^i~ decay the point of which is almost free of hadronic unccrtaintics(~ 5%). In 
what follows, we consider 200 < 1/R < lOOOGeV and analyse the dependency of branch- 
ing ratios in terms of inverse of the compactification radius R. Also, physical observables 
are sensitive to rhc where we use two different values of rhc — 0.22 and rhc — 0.29 in our 
numerical calculations. 

Figs. (l)-(4) depict the dependence of the total branching ratio in terms of the compat- 
ification parameter in two different values of rhc = 0.22 and tjIc = 0.29 for B — > TTi^i~ 
and B pi'^i~ decays, respectively. Looking at these figures, we can see: 

• Br strongly depends on the compactification radius R for both p and r channels. 
Furthermore, Br is decreasing function of the 1/R where at large value of the 
result of SM and ACD is almost the same as it is expected. Moreover, while the 
branching ratio of B ^ ni~^i~ decay is sensitive to the value of rhc for At lepton, the 
B pi^i~ decay depicts a weak dependency on the value of rhc for both p and r 
leptons. 

To sum up, we presented the analysis of the B 'ri{p)fr£^ decay in ACD model with 
single universal extra dimension. The only free parameter of the model is compactification 
radius 1/R. We studied the dependence of the branching ratio on the inverse of com- 
pactification radius 1/R. The value of the branching ratio was obtained larger than the 
corresponding SM value. 




(37) 
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Figure captions 

Fig. (1) The dependence of the total branching ratio oi B ^ tt/x+zx" on 1/R for 
rhc = 0.22, 0.29. 

Fig. (2) The same as in Fig. (1), but for the r lepton. 

Fig. (3) The dependence of the total branching ratio of S — > pii^ii~ on 1/R for rhc — 
0.22, 0.29. 

Fig. (4) The same as in Fig. (3), but for the r lepton. 
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